Spin-liquid versus dimerized ground states in a frustrated Heisenberg antiferromagnet 
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We present a Density Matrix Renormalization Group (DMRG) study of the ground-state proper- 
ties of spin-1/2 frustrated Ji — J3 Heisenberg n;-leg ladders (with ni up to 8). For strong frustration 
(J3/J1 — 0.5), both even- and odd- leg ladders display a finite gap to spin excitations, which we 
argue remains finite in the two-dimensional limit. In this regime, on odd- leg ladders the ground 
state is spontaneously dimerized, in agreement with the Lieb-Schultz-Mattis prediction, while on 
even-leg ladders the dimer correlations decay exponentially. The magnitude of the dimer order 
parameter decreases as the number of legs increases, consistent with a two-dimensional spin-liquid 
ground state. 
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Despite many years of intense investigations, the ex- 
istence of a homogeneous spin-liquid ground state for 
a spin-1/2 system on a two-dimensional square lattice 
remains controversial. This is mainly because there 
has been no definite evidence so far that a microscopic 
model could stabilize a homogeneous non-magnetic phase 
with one electron per unit cell. In fact, the known 
non-magnetic phases of spin-1/2 quantum antiferromag- 
nets in one or two dimensions usually display a spon- 
taneously broken translation symmetry related to spin- 
Peierls dimerization, as in the frustrated Heisenberg 
chain [HQ and ring exchange models, 0, or arising from 
a doubling of the ground-state unit cell, as in the Heisen- 
berg two- leg ladder. 

For one-dimensional systems a rigorous result, the 
Lieb-Schultz-Mattis (LSM) theorem, [j| implies that a 
gapped non-magnetic phase is in general associated with 
a broken translation symmetry. This result can be also 
extended to spin-1/2 Heisenberg models defined on odd- 
leg ladder geometries. There have been several re- 
cent attempts to generalize the LSM result to two di- 
mensions 0, S 01- Here it has been argued that the 
gapped phase is associated with a ground-state degen- 
eracy. However, there are different opinions on whether 
this degeneracy necessarily implies a spontaneously bro- 
ken translation symmetry (and thus the non-existence 
of a two-dimensional spin-liquid) or whether it is as- 
sociated with a topological degeneracy of fractionalized 
spin-liquid phases p], llil llll LL2| . Recentl y, o n the basis 
of a variational approach, Sorclla et a/.,[13j have pro- 
posed that a spin-liquid ground state can be stabilized 
in two-dimensions and yet satisfy the constraint imposed 
by the LSM theorem. In fact, within the formalism of 
projected BCS wave functions it is possible to construct 
a gapped state which displays spontaneous dimerization 
on any odd-leg ladder thus satisfying the LSM theorem, 
but with no dimerization for even- leg ladders. The two- 
dimensional thermodynamic limit is consistently reached 
for a large number of legs since the dimer order parame- 
ter on odd-leg ladders vanishes in this limit, thus leading 



to a homogeneous spin liquid. In this case, therefore, 
the ground-state degeneracy predicted by the general- 
izations of the LSM theorem Q is not connected to a 
spontaneously broken translation symmetry but rather 
to a topological degeneracy of fractionalized resonating 
valence bond states. [12], llil 

In this paper, we examine n;-leg frustrated Heisenberg 
ladders with Hamiltonian 
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Here Si are spin-1/2 operators on a square lattice, and 
J r , J 3 are the nearest- and third- nearest neighbor an- 
tiferromagnetic couplings along the two coordinate axes. 
In the following, we use the numerical density matrix 
renormalization group (DMRG) 14| to study the ground 
state of this Hamiltonian on ladder systems with m legs 
of length L with open boundary conditions. In our cal- 
culations, we typically performed 15-20 sweeps of the lat- 
tice, keeping a maximum of m ~ 2000 states and obtain- 
ing discarded weights smaller than ~5x 10~ 7 . Our plan 
is to carry out DMRG calculations for ladders with dif- 
ferent number ni of legs. Then, by extrapolating in the 
length L of the ladders and looking at the behavior of the 
odd-and even leg systems for modest value of n; we seek 
to gain insight into the behavior of the two-dimensional 
system. 

The classical ground state of the J1 — J3 Hamiltonian 
in two dimensions displays conventional Neel order for 
J3/J1 < 0.25. For J3/J1 > 0.25 the ground state 
has incommensurate antiferromagnetic order with a pitch 
vector depending on the frustration ratio, assuming the 
value Q = (27r/3, 27r/3) at J3/J1 = 0.5, and approaching 
Q = (7r/2,7r/2), corresponding to four decoupled Neel 
lattices, for J3/J1 — > 00. For the case of quantum spin- 
1/2, in two-dimensions, the ground state is expected to 
display long-range Neel order for J3/J1 — > 0, and nu- 
merical calculations on lattices up to 32 sites suggest 
that a non-magnetic ground state could be stabilized in 
the regime of strong frustration J3/J1 ~ 0.5. (l5| This 
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FIG. 1: Spin gap as a function of the ratio J3/J1. Left: 3-leg 
ladder for N = 12 (empty triangles), 16 (full triangles), and 
20 (stars). Right: 4-leg ladder for N = 10 (empty squares), 
12 (full squares), and 16 (stars). 



work also found signatures of dimerization for values of 
J3/J1 — 0.7 in agreement with the predictions of series 
expansions. 0| 

The effects of frustration on the antiferromagnctic cor- 
relations of the ground state can be investigated by 
studying the behavior of the spin-gap as a function of 
J3/J1. In fact, in a system with long-range spin corre- 
lations, spin excitations arc necessarily gaplcss. Instead, 
if the energy cost of the lowest triplet excitation, A, re- 
mains finite in the thermodynamic limit the ground state 
has short-range spin-correlations. As shown in Fig.^ for 
ni = 3 and 4, the spin gap increases as the frustration 
ratio J3/J1 increases. This is seen for both the even- and 
the odd-leg ladder case. In particular, in the odd-leg case, 
which is gaplcss with pow er-law spin correlations in the 
pure Heiscnbcrg limit. 1 171] the spin gap due to the finite 
length of the ladder remains almost constant for small 
values of J3/J1 but increases sharply for J3/J1 ~ 0.4, 
suggesting a transition to a gapped non-magnetic phase. 
Alternatively, in the even-leg ladder case, a finite correla- 
tion length is expected for small J3/J1 and the spin-gap 
increases smoothly with J3/J1 as no magnetic transition 
is expected. In both cases, the spin gap reaches a max- 
imum for intermediate values of J3/J1 where the effects 
of frustration are expected to be the strongest, then it 
decreases again for large J3/J1 when the limit of four 
decoupled Heiscnbcrg lattices is eventually recovered. 

The size scaling of the spin gap is shown in Fig. [2] 
For weak frustration, the spin gap extrapolates to zero 
for the odd-leg ladders, and to a constant, which de- 
creases with ni, for the even-leg ladders. [Tsj This is con- 
sistent with the gapless Neel ordered phase expected in 
the two-dimensional limit. Instead, for J3/J1 = 0.5 the 
spin gap extrapolates to a constant as L — > 00 for both 
the even- and the odd-ladders we have studied. The dif- 
ference between the regime of low and high frustration 
is also seen from the dependence of the spin-gap on the 
number of legs for a fixed chain length L (see also Fig. [2J- 
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FIG. 2: Size scaling of the spin gap for J3/J1 = 0.1 [(a) and 
(b)], and J3/J1 = 0.5 [(c) and (d)]. (a) and (c): spin gap as a 
function of the length of the ladders, L, for different number 
of legs, ni. (b) and (d): spin gap as a function of m, for 
L — 10. Empty (full) symbols correspond to odd- (even-) leg 
ladders. Lines are guides for the eye. 



(b) and (d)). For low- frustration the spin-gap decreases 
with the number of legs both for even (full symbols) and 
odd (empty symbols) n;. However, in the regime of high 
frustration it decreases with n/ only for even leg samples 
while it increases with the number of legs on odd-leg sam- 
ples. This behavior is consistent with a two-dimensional 
phase which has exponentially decaying spin correlations. 

The presence of a finite-gap in the excitation spectrum 
of the nj-leg ladders has consequences in view of the LSM 
theorem. In fact, on odd- leg ladder systems it is possi- 
ble to construct an excitation in the singlet sector with 
momentum (jt, 0) which becomes degenerate with the 
ground state in the thermodynamic limit. This implies 
cither a gaplcss spectrum or, in the presence of a finite 
gap, a two-fold degenerate ground state with a doubling 
of the unit cell and a spontaneously broken translation 
symmetry In the one dimensional model this is known to 
be realized through spin-Peierls dimerization. Instead, 
the LSM result does not apply to even-leg ladders so 
that, in these geometries, both translationally invariant 
and dimerized ground states are in principle compatible 
with a finite triplet gap. 

The occurrence of spin-Peierls dimerization can be 
studied by calculating the response of the system to 
a nearest-neighbor spin-spin operator which breaks the 
translation symmetry along the chains with momentum 

Q = (7T,0) 

d = J2e iQr S r -S r+x ■ (2) 

r 

Here x = (1, 0) is a unit vector along the chain direction. 
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FIG. 3: Upper panels: D(S) for J3/J1 = 0.5. Left: one 
dimensional (Majumdar-Gosh) chain for N — 12 (triangles), 
and 24 (squares). Right: 2-leg ladder for N = 12 (triangles), 
and 30 (squares). Full (empty) symbols correspond to 5 < 
(S > 0). Lower panel: size-scaling of the dimer order param- 
eter, D(S — > 0), as a function of the length of the ladders, L, 
for different numbers of legs, n;. 



This can be done in general by adding to the Hamilto- 
nian a term SO where 5 is a (small) generalized field. The 
order parameter can be calculated as the limit for d — > 
of the ground-state expectation value of O in presence of 
the field, D = ]ims^>o{0)s/N, where N = L x ni is the 
total number of sites. On periodic finite-size systems D 
vanishes in general by symmetry as it breaks the trans- 
lational invariance and the symmetry under site- centered 
lattice reflections along the chain direction of the unper- 
turbed Hamiltonian Q. However, on samples with an 
even length L and open boundary conditions, the dimer 
order parameter (0) will be in general non-zero and can 
be calculated using the Hellmann-Feynman theorem as 
D = de(8)/d5\s=o- Here e(S) is the ground-state energy 
per site (in unit of Ji) in the presence of the perturbation. 
As a result, within the DMRG technique the dimer order 
parameter can be calculated with simple energy measure- 
ments by computing e(<5) for a few values of <5 and then es- 
timating numerically the limit D = lim,5^o(e(<5) — eo)/S. 
This is illustrated in the upper panels of Fig.|3|for a single 
chain and a two-leg ladder at J3/J1 = 0.5. Here, as a con- 
sistency check, the dimer order parameter is estimated by 
calculating the limit for 5 — > of D(6) = (e(S) — eo)/S 
both for positive (filled symbols) and negative (empty 
symbols) <5's. The two limits converging to the same 
value. In particular, for the one-dimensional chain at 
the exactly solvable point J3/J1 = 0.5 (Majumdar-Gosh 
model) the known size-independent result, D — 0.375, is 
recovered Q. 

The size scaling of the dimer order parameter obtained 
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FIG. 4: Dimer susceptibilities for the 2-leg (left panels) and 
3-leg (right panels) ladders. The lower panels show x($)/N 
vs S for various lengths, and the top panels the size scaling 
of the extrapolated values, % — lim«_o x($) ( see text). Note 
the different normalizations for the 2- and 3-leg cases in the 
upper panels. Lines are guides for the eye. 



with this procedure is shown in the same figure for lad- 
ders with different numbers of legs. This analysis reveals 
close similarities with the variational scenario of Ref. . 
In fact, in the odd- leg ladder case the dimer order param- 
eter extrapolates to a constant for infinite chain length, 
as required by the LSM theorem. Instead, in the even-leg 
cases, where the system is not constrained by the LSM 
theorem to dimerize, the dimer order parameter extrap- 
olates to zero. 

This conclusion is supported also by the calculation 
of dimer susceptibilities. These can be calculated within 
our numerical approach by considering ladders with odd 
length L. In these geometries the unperturbed Hamil- 
tonian is reflection symmetric around the central col- 
umn of sites, so that the order parameter D vanishes 
by symmetry for any finite-size cluster, and the ground- 
state energy has corrections proportional to S 2 . There- 
fore e(5) ~ eo — x<5 2 /2 , with \ tne generalized suscep- 
tibility associated with the operator O, namely, x = 
2(il) \6{E [) - H)- 1 d\ip )/N. If true long-range order 
in the dimer correlations exists in the thermodynamic 
ground state, the finite-size susceptibility will diverge as 
the system size increases. In particular, it can be shown 
that it is bounded from below by the system volume 
squared, x ~ N 2 . 0] Thus susceptibilities are a sensi- 
tive tool for detecting the occurrence of long-range order. 
In analogy with the calculation of the order parameter, 
the susceptibility x = —d 2 e(S)/dd 2 \s=o can be calculated 
numerically from x = lim^o x(<$) = — 2(e(5) — eo)/S 2 , 
as illustrated in the bottom panels of Fig. 0] for 2- and 
3-lcg ladders. 
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FIG. 5: Left: Size scaling of x/N 2 for the 3- and 5-leg lad- 
ders as a function of the chain-length L. Right bottom: Size 
scaling of \ f° r the 2-, 4- and 6-leg ladders as a function of 
the chain-length L (bottom). Right top: Size scaling of the 
L — > oo extrapolated values of the even-leg ladder suscepti- 
bilities, Xoo, as a function of the numbers of legs, n;. 



The behavior of the susceptibilities for the even- and 
odd-leg ladders is remarkably different (Fig.^J. In partic- 
ular, x/N decreases with the linear size L in the even-leg 
ladder case, and it increases with L in the odd-leg ladder 
case. The susceptibilities for the odd-leg ladders appear 
to diverge as N 2 , as required for true long-range dimer 
order. In contrast, for the even-leg ladders the suscepti- 
bilities are bounded, indicating short-range dimer corre- 
lations with a finite correlation length. 

In order for the two-dimensional limit to exist, the 
ground-state correlations for even and odd-leg ladders 
must converge in the limit of a large number of legs. 
This is the case, for instance, in the limit of small frus- 
tration of this model, where the spin correlations decay 
exponentially for even ni and with a power-law for odd 
ni . Here the two-dimensional limit is reached as the cor- 
relation length for the even-leg ladders diverges expo- 
nentially with ni leading to long-range antifcrromagnctic 
order in two dimensions. [TtI llq As we have shown, in 
the regime of strong frustration, the dimer-corrclations 
are short-ranged on even-leg ladders while a finite-dimcr 
order parameter is observed on odd- length ladders. How- 
ever, as the number of chains ni is increased, the odd-leg 
ladder dimer order parameter decreases (Fig. [3J , and the 
divergence of the dimer susceptibility becomes weaker 
(see left panel of Fig.[SJ). On the other hand, the infmite- 
L dimer susceptibility, Xooi on even- leg ladders does not 
appear to diverge as the square of the number of chains, 



lation symmetry in the two-dimensional limit appears as 
a plausible interpretation of our results. 

In conclusion, we have shown that a spin-gapped 
ground state with short-range antiferromagnetic corre- 
lations is stabilized by frustration on the spin-1/2 Jy — J3 
model at J3/J1 = 0.5 on ladders with n; = 1 to 8 legs. 
The behavior of the spin gap by increasing the number of 
legs is consistent with a non-magnetic ground state in the 
two-dimensional limit. On odd- leg ladders we find a finite 
dimer order parameter associated with a spontaneously 
broken translation symmetry. However, as the number of 
legs increases n; = 1, 3, 5, 7 the size of the order param- 
eter decreases and the divergence of the associated sus- 
ceptibility becomes weaker. These results suggests that 
in the two-dimensional limit the dimer order parameter 
vanishes. Although the numerical data we have presented 
here was for J3/ Ji = 0.5, we find similar results for other 
values of J3/J1 near 0.5 and believe that for a range of 
J3/J1 values this model exhibits a non-magnetic ground 
state in two dimensions. Our results are consistent with a 
recently proposed scenario where the odd-leg dimer 
order vanishes for infinite number of legs leading to a 
homogeneous spin-liquid in two dimensions. 
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(see right-panels of Fig. [SJ as one would expect in 
presence of long-range dimer order in two dimensions. 
Thus, a ground-state with no spontaneous broken trans- 
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